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^ ' Abstract 

The spectra of octet and decuplet baryons are studied using SU{'i) lattice QCD at the 
' quenched level. As an implementation to reduce the statistical fluctuation, we employ the 

JnI , anisotropic lattice with 0{a) improved quark action. In relation to A(I405), we measure 

also the mass of the SU{?)) flavor-singlet negative-parity baryon, which is described as a 
^— ^ , three quark state in the quenched lattice QCD, and its lowest mass is measured about 1.6 

GeV. Since the experimentally observed negative-parity baryon A(1405) is much lighter 
, than 1.6 GeV, A(1405) may include a large component of a NK bound state rather than 

the three quark state. The mass splitting between the octet and the decuplet baryons are 
Oh, also discussed in terms of the current quark mass. 

CD 

1 Introduction 

;_j ■ The lattice QCD simulation has become a powerful method to investigate hadron properties 

directly based on QCD. The hadron spectroscopy in the quenched approximation, i.e. without 
the dynamical quark effect, has been almost established, and reproduces experimental values of 
the low-lying hadron masses within 10 % deviation Q. Detailed investigation with dynamical 
quarks will give us an insight on the dynamical quark effect on hadron properties |^]. With 
these situations, we now have a stage to proceed into more extensive and systematic studies of 
hadron structure in lattice QCD, and compare the lattice results with the model analysis or the 
phenomenological approaches such as the potential models The latter point of view would 
give us clearer physical picture of extracted information from lattice simulations, and may be 
useful to extract the physical quantities beyond the lattice QCD applicability. Furthermore, 
comparing hadron properties in the lattice QCD with the potential model analysis using the 
potential derived from lattice QCD, one can verify the applicability of these approaches in a 
self-consistent manner. This is our motivation of lattice calculation of the static three quark 
potential, which is responsible to the baryon properties Q . 

Compared with the ground state hadrons, excited state hadrons are far from established. 
The purpose of our present study is to perform detailed investigation of excited states as well 
as the ground state hadrons. Among them, in this paper, we focus on two important subjects: 
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(a) Origin of Octet-Decuplet baryon mass splitting. 

There have been proposed several models to explain the low-lying hadron spectrum. As 
one of celebrated models, the nonrelativistic quark model explains the N-A mass splitting 
with one-gluon-exchange (OGE) interaction and results in ^ 

where Mj is the constituent quark mass. This implies mass difference between octet and 
decuplet baryons decreases with increasing quark mass. In lattice QCD simulation, it is 
possible to change the current quark mass through the hopping parameter k. Then we 
compare the obtained mass splitting with the form (1), although the relation between 
current and constituent quark masses is not so clear. 

(b) Structure of A(1405) and other negative parity baryons. 

The detailed lattice study of the negative-parity baryons has been started rather recently 
0. Among the negative-parity baryons, we pay much attention to A(1405) with = 
l/2~, since its structure in terms of the constituent quark picture is not well understood. 
There are interesting two possibilities proposed for A(1405), an SU{3)f singlet state 
{qqq) and a NK bound state as (qqq-qq)- In this paper, we investigate the flavor-singlet 
baryon spectrum with spin 1/2 and negative parity in lattice QCD, and compare it 
with A(1405). The second possibility as the NK bound state is now in progress at the 
quenched level, where quark-antiquark pair creation (strictly speaking, dynamical quark 
loop effect) is absent and then the quark-level constitution in hadrons is definitely clear 
in the simulation. 

In lattice QCD simulations, not only the quark mass, we can also change the number 
of dynamical quark flavor. This enables us to extract the quark loop effect on the hadron 
properties. However, the flavor-number dependence does not seem significant for the low-lying 
hadron spectrum. Then, we start with calculations at the quenched level (with no dynamical 
flavor). Before proceeding to the numerical calculation, however, we should describe a technical 
problem and equipment to circumvent it. 

The difficulty of extracting the excited state masses lies in the rapid growth of the statistical 
fluctuations in the correlation functions. In the practical simulation for the correlation function, 
it is hard to identify the reliable range of its Euclidean temporal distance where the relevant 
information is kept without suffering from the large statistical noise. To overcome this problem, 
we adopt the anisotropic lattice, on which the temporal lattice spacing Ur is finer than the 
spatial one, Ua 0. Detailed information in the temporal direction makes these analyses as the 
mass measurements extensively easy. This approach is especially efficient for the heavy particle 
correlators, such as excited states and glueballs, for which the noises grow rapidly against the 
signals. On the other hand, introduction of anisotropy requires us additional effort in the 
numerical simulations. Due to the quantum effect, the renormalized anisotropy ^ = aa/ar 
differs from the bare one in general, and at the first stage of the simulation one need to tune 
the bare anisotropy so that the quark field retains the same renormalized anisotropy as the 
gluon field. 

This paper is organized as follows. In the next section, we briefiy summarize the anisotropic 
lattice, especially the quark action. Then, we describe the numerical simulation and discuss on 
the obtained result. The last section gives our summary and outlook. 
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2 Anisotropic lattice 



The anisotropic lattice has become an extensively useful tool in the lattice QCD simulations. 
In addition to aforementioned advantage, fine temporal resolution is particularly significant to 
extract the information of hadron correlators at finite temperature. Another advantage is that 
it can treat the relatively heavy quark without introducing the effective theoretical approaches. 
This feature is particularly suited for the study of charmonium and charmed hadrons. Although 
this is not a subject of this paper, this is one of reasons we adopt the anisotropic lattice. In 
this work, we treat the hadron correlators, and hence we focus here on the quark action on the 
anisotropic lattice. 

As the quark action on the anisotropic lattice, we adopt the 0{a) improved Wilson action 
1^]. The construction of the action is along the program of Fermilab formulation Their 
treatment incorporates the full quark-mass dependence in the construction of lattice quark 
action, so that the quark mass region with ~ is also available without large 0{ma) 
uncertainty. The Fermilab approach is naturally generalized to the anisotropic lattice. There 
is some arbitrariness in choosing the parameterization of action, and we use the form proposed 
in 12|. In |12|, the advantages of this form are discussed in detail. The quark action is 
written in the following form: 

SF = Y.^{x)K{x,y)i;{y), (2) 

K{x,y) = 6,,y - Kr{{l-j4)U4{x)6^+iy{l + j4)ul{x-4)6^_^^y} 

- ^^T.{(^-^i)U^{x)S,+ly{r + l^)U|{x-^)5^_iJ (3) 

i 

i ij 

The gluon field is represented with the link variable ~ exp(—iga^Af^), and ip denotes the 
anticommuting quark field. The spatial and the temporal hopping parameters. Kg- and k^, 
respectively, are related to the bare quark mass mg and the bare anisotropy parameter as 

Ka- = l/2(mo +7F + 3r), = 'JpHa, (4) 

where mo is measured in the spatial lattice unit. The value of the Wilson parameter r is set as 
r = 1/^. The coefficients ce and cb in the clover terms are introduced to eliminate the 0{a) 
error induced by the Wilson term, and coincide with unity at the tree level. 

We apply the mean-field improvement which reduces large contributions from the tadpole 
diagrams to the renormalization. This is achieved by replacing the link variable as [/ — > U/uq, 
with the mean-field value of the link variable, uq. On the anisotropic lattice, two mean- field 
values Ua and Ur are defined for the spatial and the temporal link variables, respectively. We 
employ the definition of the mean-field value through the average of U in the Landau gauge. 
With Ua and Ut, the mean- field improved values of the clover coefficients at the tree level are 
expressed as c^; = l/n^Uo- and cb = ^/u^- The anisotropy parameter is related to the improved 
anisotropy as jf = If • u^^/ur- It is convenient to define k, as 

- = ^-2(7F + 3r-4) (=2(mo + 4)). (5) 

For the light quark systems, the extrapolation to the chiral limit is performed in 1/k. 
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In practical simulations, the anisotropy parameter should be tuned so that the fermionic 
anisotropy defined with fermionic observable coincides with the gauge field anisotropy. This 
is called as "calibration" . Several procedures have been used for the calibration. In this paper, 
we set the value of 7^ using the dispersion relation of the pseudoscalar and the vector mesons. 
We assume that the meson field is described with the lattice Klein-Gordon equation, 

S = ^Y.<t^''^^)\-^FDl-& + mt\m, (6) 

with the lattice covariant derivative. The free meson field satisfies the dispersion relation, 

cosh. E{p) - cosh£;(0) = f 12^%. (7) 
This relation is used to define the fermionic anisotropy W^ - 

3 Numerical Results 

Lattice Setup The SU(3) lattice QCD simulations are performed on an anisotropic lattice 
of the size 12^ x 96 with anisotropy = 4, at the quenched level. As the gauge field action, we 
adopt the anisotropic Wilson action with the parameters 7^) = (5.75,3.072), determined 



by Klassen so as to give the renormalized anisotropy ^ = 4 within 1 % uncertainty [13|. At 
these parameters, the lattice cutoff defined by setting the string tension ^/a to be 427 MeV 
is found to be a~^{^/a) ~ 1.0 GeV. The gauge configuration is fixed to the Coulomb gauge, 
which is convenient to smear the quark propagators by extending the quark source spatially 
on a time slice. 

The mean-field values are determined on the lattice of half size in the temporal direction, 
12^ X 48, at the same (/3, 7g)- We fix these gauge configurations to the Landau gauge and 



determine the mean- field values and Ur self-consistently as described in |11|. They result in 
= 0.7620(2) and Uo- = 0.9871 (error is less than the last digit). 

Calibration The quark field calibration is performed along the course described in the last 
section. The pseudoscalar and the vector meson correlators are calculated with momentum 
p = 0, 27r/16 and 2 • 27r/16. We use the meson operators listed in Table § and a standard 
procedure to extract the meson energy. 

Figure |^ shows the result of the calibration. The left three values of k correspond to 
the quark masses around the strange quark mass. We use these three values of k for the 
following analysis of hadron spectroscopy. As is clearly observed in Fig. |l|, in the light quark 
region {niq ~ nis), one can set the bare anisotropy 7^? = 4. In Table ||, we list the values of 
K, the pseudoscalar and the vector meson masses for the degenerate quark case. The chiral 
extrapolation is carried out linearly in 1/k, and results in the critical hopping parameter as 
= 0.12637(2). 

Baryon spectrum As listed in Table ^, we use the standard baryon operators which have 
the same quantum numbers as the corresponding baryons and survive in the nonrelativistic 
limit. At large t (and large Nt — t), the baryon correlators are represented as 

GB{t) = Y.(B{x,t)B{x,0)) = (1+74) [cb+ ■ e-*™^+ + bcj,- ■ e-^^'-*)"^^- 

X 

+ (1 - 74) \bcB+ ■ e-(^*-*)™s+ + cb- ■ e-'^B- 1 , (8) 
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Figure 1: Results of the calibration in an anisotropic lattice. The value of at which = C 
holds is determined for each value of n. 



K rnps my 
0.124 0.1504( 6) 0.2240(22) 
0.123 0.2036( 6) 0.2602(13) 
0.122 0.2499( 5) 0.2966( 9) 



Table 1: The hopping parameters used in the spectroscopy and the PS and V meson masses 
for the degenerate quark case. 

where h = +1 and —1 for the periodic and antiperiodic temporal boundary conditions for the 
quark fields. Thus, combining the parity-projected correlators under two boundary conditions, 
one can single out the positive and negative parity baryon states with corresponding masses 
mB+ and uib-, respectively. 

In our calculation, u, d current quark masses are taken to be the same value as = rud 
or Ku = Kd, and the strange current quark mass or Ks is taken to be an independent value. 
Then, the baryon masses are expressed as the function of m„ and rUg like MB{mu,ms), and 
therefore the baryon masses MB{mu,ms) are to be depicted on the (m„,ms) plane. However, 
the lattice QCD result for the baryon masses Mb seem to be well described as a function of 
the averaged current quark mass over the three quarks, niq = (mi + m2 + 7713)/ 3, like MB{mq) 



Meson Pseudoscalax M{K) = s"f;,u 

Vector Mk{K*) = s-fku 

Baryon Octet i3Q,(S'^) ~ fjj[ua{dfjSj — spd^) — daispu^ — u/ss^)] 

Octet (A) -Bq:(A) = (0^5) i3-y[ua{di3Sj — Sfjd-y) + da{sf}U-y — ups-y) — 2sa{upd^ — d/ju-y)] 
Singlet -Bo(Ai) = {C'j5)f3'y[ua{dfjSj - spd^) + daispu-f - ups^) + Sa{upd^ - dpu^)] 

Decuplet ■Bafc(I]*°) = {C^k)p'y[Ua{dpS^ + Sj3d.y) + da{Sf3U^ + Mffg-y) + Sajupdj + d/3U^)] 

Table 2: Examples of the hadron operators. For baryon operators, the contraction with the 
color index is omitted. 
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Figure 2: The spectra of octet and decuplet baryons with positive parity. The horizontal axis 
denotes the averaged current quark mass. 

in each channel. Therefore, we will use the simplified figures as the function of rn-g, although 
the actual lattice QCD calculations are performed with different quark masses on mu{= rrid) 
and rag. 

Let us start with the ground-state baryons shown in Figure |2|, where the horizontal axis 
denotes the averaged current quark mass over the three quarks in the physical unit. We use 
the naive relation between the current quark mass and the hopping parameter as 



1. 

mq = -{mi +m2 + ms) 



mi 



1 



(9) 



The baryon masses are linearly extrapolated to the chiral limit, mq = 0. In the analysis of 
the baryon spectrum, we use the lattice scale determined by setting the averaged mass over 
the octet and the decuplet baryons at the chiral limit to the averaged mass of and A. It 
results in = 0.9 GeV. The strange quark mass is determined so that the octet and the 
decuplet baryons in Fig. |2| globally reproduce the experimentally measured masses. Thus, for 
the strange quark, we adopt Kg = 0.1229(1), which roughly corresponds to ~ 0.16 GeV as 
the current quark mass. This value seems consistent with the standard strange current-quark 
mass. Hereafter, we will fix Ks = 0.1229(1) for the strange quark. Figure |2| shows that the 
lattice result reproduces that the measured baryon spectrum within 5 % deviation. 



Octet-Decuplet mass splitting Let us consider the difference of octet and decuplet baryon 
masses. This mass splitting monotonously decreases as the current quark mass increases. This 
tendency is consistent with the one-gluon-exchange explanation in the nonrelativistic quark 
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Figure 3: The positive and negative parity baryon masses in the flavor octet and singlet 
channels. 



model, assuming the constituent quark mass Mq increases with the current quark mass rriq, as 
Mg ~ Mo + rriq. 

Negative parity baryons Figure ^ shows the mass spectra of the flavor octet and singlet 
baryons with positive and negative parities as the function of the averaged current quark mass 
niq. In each channel, the baryon mass spectrum shows clear linear behavior on the averaged 
current quark mass niq in the whole measured mass region. Then, we extrapolate the masses 
to the physical situation m„ = rud — 0, = 0.16GeV, i.e., niq ~ nis/S. One striking feature 
is that the singlet and octet negative parity baryons are almost degenerate. On the other 
hand, the lowest mass of the flavor-singlet positive-parity baryon is much larger than those in 
other channels. Extrapolating to the physical situation m„ = rrid ~ 0, = 0.16GeV {niq ~ 
nis/3), we obtain 1.6 GeV for the lowest mass of the flavor-singlet negative-parity baryon, 
which is described as a three quark state in the quenched lattice. Then, the experimentally 
observed negative-parity baryon A(1405) is much lighter than the lowest mass of the flavor- 
singlet negative-parity baryon obtained in lattice QCD at the quenched level. Since the q-q 
pair creation is absent in the quenched QCD, this result implies that the simple three valence 
quark picture would not valid for A(1405). One possible explanation is, as frequently suggested, 
that the A(1405) is a mixture of a three quark state and an NK bound state. To clarify this 
subject, we are going to measure the NK state on the quenched lattice, where genuine bound 
state properties is apparent owing to the absence of the q-q pair creation. 
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4 Summary and outlook 



In this paper, we report present status of our investigation of hadron properties using lattice 
QCD simulations. At this stage, systematic investigation of baryon spectrum is in progress at 
the quenched level. The result on the anisotropic lattice with ~ 0.9 GeV and anisotropy 
^ = 4 is as follows: (a) The octet-decuplet baryon mass splitting decreases with increasing 
current quark mass. This tendency is consistent with the one-gluon-exchange explanation in 
the constituent quark model, (b) The negative-parity baryons in octet and singlet channels are 
measured in good statistical precision. However, the experimentally observed negative-parity 
baryon A(1405) seems much lighter than the lattice QCD result for the lowest mass of the 
SU(3) flavor-singlet negative-parity baryon. This may suggest the possibility of that A(1405) 
is a mixture of the NK bound state and the three quark state. This would be clarified by 
successive lattice calculations for qqq-qq system in the quenched approximation. 

Our present results have been carried out on rather course lattice. Hence, we need to 
perform the simulations on finer lattices to remove lattice artifacts. Then, the anisotropic 
lattice will be a powerful device to study detailed properties of hadrons including excited 
states, exotics and glueballs. The simulation with dynamical quarks are also to be performed. 

Another course of our program is a comparison with the potential model analysis using the 
static quark potential extracted from lattice QCD. This will give us important information on 
the quark wave function and novel insight on the quark structure of hadrons. 

We are grateful to Profs. Il-Tong Cheon and Su Houng Lee for their warm hospitality at 
Yonsei University. The lattice QCD simulations have been performed on NEC SX4 at Osaka 
University, and Hitachi SR8000 at KEK. 
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